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Abstract 

We obtain the exact solutions for a family of spin-boson systems. This is achieved 
Q>^ . through application of the representation theory for polynomial deformations of the 

su{2) Lie algebra. We demonstrate that the family of Hamiltonians includes, as 
■ special cases, known physical models which are the two-site Bose-Hubbard model, 

the Lipkin-Meshkov-Glick model, the molecular asymmetric rigid rotor, the Tavis- 
Cummings model, and a two-mode generalisation of the Tavis-Cummings model. 
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m 1 Introduction 
od 

^ . The study of polynomial deformations of Lie algebras is an area of research which has 
found many applications in systems involving non-linear interactions pT-TS]. In recent 
publications |4j we have formulated such methods for the analysis of a class of multi- 
X : boson systems. The approach of 0] is to express the Hamiltonian of the systems in 
5y , terms of the generators of polynomial deformations of the sit (2) Lie algebra, through the 
explicit construction of Fock-space representations. By utilising a correspondence between 
the Fock-space representations and differential operator realizations, it was shown that 
exact solutions are obtained in terms of a system of coupled equations. These equations 
can be viewed as providing a Bethe ansatz type of solution for the calculation of the 
energy spectrum and associated eigenstates. The generality of this approach allows for 
application on a wider level. The work described below is concerned with extending these 
methods to the study of a family of Hamiltonians which couple multi-boson degrees of 
freedom to a spin degree of freedom. In this manner we unify the problem of exactly 
solving spin-boson Hamiltonians to a particular class which contains within it a number 
of models which are already known in the literature, as we will discuss. 
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The main result of this paper is the derivation of the exact eigenfunctions and energy 
eigenvalues of the infinite family of spin-boson systems defined by the Hamiltonian 

M 

H = Y,w^N^ + g'4 + g[Jla\'■■■al'' + raf'---al'''^''^, (1.1) 
1=1 

where throughout r, s G Z+, M,ki, - ■ ■ ,kM G N, a,,, a] and Ni = a|ai are bosonic an- 
nihilation, creation, and number operators respectively, J±fi are the generators for the 
su{2) spin algebra, and Wi, Wij, g are real coupling constants. The Hamiltonian of the 
form (11. ip appears in the description of various physical systems of interest in atomic, 
molecular, nuclear and optical physics. We will explicitly demonstrate that (II .ip includes 
as special cases several known models which are the two-site Bose-Hubbard model p] , the 
Lipkin-Meshkov-Glick model [6], the molecular asymmetric rigid rotor [7j, the many-atom 
Tavis-Cummings model [8l|9] and a two-mode generalized Tavis-Cummings model [TO] . 

This paper is organized as follows. In section 2 we introduce new higher order polyno- 
mial algebras which are dynamical symmetry algebras of Hamiltonian (II. ip , which enable 
an algebraization of the spin-boson systems. We construct finite-dimensional unitary 
representations of the dynamical symmetry algebras in section 3 and the corresponding 
single-variable differential operator realizations in section 4. This leads to the higher or- 
der differential operator realizations of the Hamiltonian (II. ip . In section 5 we establish 
the quasi-exact solvability of this differential operator [TlHTS] and solve for the eigenvalue 
problem via the functional Bethe ansatz method (see e.g. [HHIE]). In section 6 we present 
explicit results for several special cases, thus providing a unified derivation of exact solu- 
tions to the widely-studied models mentioned above. We summarize our results in section 
7 and discuss further avenues for investigation. 



2 Algebraization 

In this section we introduce new higher order polynomial deformations of s/(2) and give 
an algebraization of the Hamiltonian (II. ip . Our approach extends previous studies [iHl] 
where more restricted classes of systems have been exactly solved using polynomial algebra 
structures. 

We introduce generators 

M M (i) 



n = nn^^^ 'p- = p-Uq?^ ^0= ° ^;=;^° , (2.1) 



i=l 1=1 



where 



Po = y, P+ = Jl, P- = J'., 
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are M + 1 mutually commuting operators. It can be shown that V±^o satisfy the following 
commutation relations: 



M 



i=l 



M 



(2.2) 



i=l 



where C is the Casimir operator of su{2), 

^ = ^^^^^^^^ > A = Q«-Qr\ ^ = l■■■M-l, 
are M central elements of (12. 2p and 



(2.3) 



ij^^'^JC,Vo,C) = -Y[[C -{rlC + rVo + r-i + l){rlC + rVo + r~i)] 



1=1 

hi 



i=l 



11=1 



are polynomial functions of degree 2r and fcj respectively. Thus (12. 2 p defines a polynomial 
algebra of degree 2r + ^^f-,^fcj — 1. 

In terms of the generators of the polynomial algebra (12. 2p . the Hamiltonian (II. ip can 
be written as 



M 



i=l 



M 



n 



i=l 



(V+ + V. 



(2.4) 



with the number operators having the following expression in ^ and £i 

( K \ \ 1 



It follows that the polynomial algebra (12. 2p is the dynamical symmetry algebra of the 
Hamiltonian (II. ip . 



3 Unitary irreducible representations 

Irreducible representations of the polynomial algebra (12. 2p can be constructed in the 
tensor product space of the representation space of P± o and the Fock spaces of {Q±''o}- 
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As shown in [1], the Fock states for irreducible representations of {Q^^^} are labelled by 
quantum numbers Qi = ^W'^ ' ' ' ' ''^'"fc^^'^^ ? through 

i i i 

^\ki (mi+qi-k~^) 

\qi,mi) = — t^?ttt|Q)- = O'l'---- (3-1) 

[ki{mi + Qi - )J! 

The action of Q±o on these states is 

QQ^\qi,mi) = [qi + mi)\qi,'mi), 



Q^l^\qi,mi) = Y[ (^^i + Qi + ^^'j^2 ^ ^ |gi,mi + l), 
Q-\qi,mi) = [[yuii + qi p j \qi,mi-l}. 



The irreducible representations of P±fi can be deduced from the S'u(2)-module Vj, j = 
0, |, 1, ■ ■ ■ as follows. First, it can be shown that Po,± satisfy the relations 

= i;^'^\Po,C)-i;^'^\Po-l,C), (3.2) 

where 

r 

^(2^) (Po, C) = -l[{C-{rPo + r-i + l)(rPo + r - «)) (3.3) 
1=1 

is a polynomial in Pq and C of degree 2r. Thus fl3.2p is a polynomial algebra of degree 
2r — 1. The Casimir operator of (13.21) takes fixed value HLil^ ~ ^(^ ~ ^))- 
It is easily verified that there are min{r, 2j + 1} lowest weight states, 

\J,0;P) ~ Jl\j,0), p = 0,l,---,min{r-l,2j}, 

where | j, 0) is the lowest weight state of su{2). This implies that finite-dimensional irre- 
ducible representations of (13. 2p . denoted as Vj^p, are labelled by quantum numbers j and 
p, j = 0, |, 1, ■ ■ ■, and p = 0,1, ■ ■ ■ , min{r — 1, 2j}. Thus we have the branching rule from 
su{2) representation Vj into Vj^p of (13. 2p : 

Vj — typ=o Vj^p. 

General basis vectors in the irreducible representation space Vj^p are given by \j, n; p) ~ 
(P+)"|j,0;p). Explicitly, 



The action of Po,± on these vectors is given by 
Po\j,n;p) = (- — - + n] \j,n;p), 



P+\j,n;p) = J]^ a/(p + 2 + rn){2j - p - i + 1 - rn) \j,n + l;p), 
1=1 

r 

P^\j,n;p) = Y\ V (p ~ ^ + 1 + rn){2j — p + i — rn) \j, n — 1; p). (3.5) 

i=l 

It can also be shown that 

P_|j,0;p) = 0, P+|j, = 0, 

r 

where A is a non negative integer taking specific values A = 0, 1, ■ ■ ■ , min{r — l,2j} 
according to j and p. Moreover, ^■'"^'^ is always a non-negative integer. Therefore n = 
0, 1, ■ ■ • , and (13. 5p is a finite-dimensional representation of (13. 2 p with dimension 

2j-p-\ I ^ 

r 

We now construct irreducible representation of (12. 2p in the tensor space V^-,p®'Hqi'' ■■■(>?) 
T-t[^\ where V,-,p is the representation space of P±,o and HqJ is the Fock space of 
From (12. 3 p we have 

M-l M-1 

Q« = Q^^ + ^ C„ (M + 1)/C = MPo + MgJ"') + 5^ fiC^. 
This implies that for any irreducible representation of (12. 2p defined by basis states |j, n; p)® 



M-l 

mi = tum + Qm - Qi + ^ IfM, i = 1, ■ ■ ■ , M - 1, 



11=1 



M-l 



where k and denote the eigenvalues of central elements JC and C^, respectively. It 



follows that qi < qm + Zlul/ ^^"^ 



. M + 1 1 \ - , v-^ , p-j 

n + mi = Ai, Ai = k - v7 2^ /^V + V " 



Clearly Ai always take non-negative integer values, i.e. = 0, 1, . . .. Thus, the irreducible 
representation of (12. 2p has basis states 

M 

\hn,p,{q},{l},n) = ® JJ 



i=l 



n (3.6) 



^{p + rn)\{2j -p-rn)\ fj^ ^[k,{Ai + q, - kf - n)\\ 
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where 

' 0,l,---,min{AM,^^} forM>0 



n 



1... "M^p^ for M = 

' r 

The action of fl2.ip on these states is given by 

f P — 3 \ 
'PQ\j,n,p,{q},{l},K) = \-n- k] \j,n,p,{q},{l},K), 



'P+\j,n,p,{q},{l},K) = JJ V(p + z + rn)(2j - p - i + 1 - rn) 

x\j,n + l,p,{q},{l},K), 



'P-\j,n,p,{q},{l},K) = YlViP - ^ + ^ + - P + ^ - 

i=l 

M kj ^ ,,h ^ \ 1/2 

X 

i=l ^=1 



nn y^i + 1i + ^2 "J \j,n-l,p,{q},{l},K). 



This gives an TV" + 1 dimensional representation of the polynomial algebra ( I2.2p . where 



min{AM, ^^^} for M > 
^teA for M = 



4 Differential operator realization 

The finite-dimensional irreducible representations in the proceeding section can be realized 
by differential operators acting on A/" + 1-dimensional space of monomials with basis 
|l, 2;, 2;^, ^-'^j, by mapping the basis vectors (I3.6P into monomials in z: 

\3^n,p,{q},{l},n) — 



\p + rn)\{2j-p- rn)\ ]\-U[h{A, + Qi - ^ - n)]\ 

The corresponding single- variable differential operator realization of (12.11) in the monomial 
space takes the following form 

-n ^ , P-j 
Vo = z- /tH , 

az r 



V+ = zY\_('2j—p — i + l—rz^j 
1=1 ^ 



i=l v=l 



ttW (rz— +p-i + l\ . (4.1) 



Note that V- contains no singularities as 111=1 (P — + 1) = for all allowed p values. 

We can thus equivalently represent Hamiltonian (12. 4p as the single- variable differential 
operator of order M. = max{r + ^i; s}, 

H = ^WiNi + g' (rz-^- j +p\ + 5^2;"^ JJ f + p - z + 1 j 

i=l ^ ' i=l ^ ' 

+9Z f 2j - p - i + 1 - rz-^j 
1=1 ^ ^ 



flf[,,U,,,-t^^-4] (4.2, 



X 

i=l u=l 



with 



M-l ^ M-1 



d M + 1 p-j I \ 1 

5 Exact solutions 

We will now solve for the Hamiltonian equation 

Hi){z) = Ei){z) (5.1) 

for the differential operator realizations by using the functional Bethe ansatz method 
[T314T6] . where iIj{z) is the eigenfunction and E is the corresponding eigenvalue. It is 
straightforward to verify 

Hz- = z-+'gf[i2j-p-z + l-rn)f[f[kJA + q,-n-^^^^^^^^) 

i=l i=lu=l ^ l^i / 

+ lower order terms, n G Z+. (5.2) 

This means that the differential operator (14. 2 p is not exactly solvable. However, it is 
quasi exactly solvable, since when n = M the first term (~ z"^^"^) on the r.h.s. of (15. 2p is 
vanishing. That is H preserves an invariant polynomial subspace of degree M , 

HVCV, V = span{l,2,...,z^} (5.3) 

Thus up to an overall factor, the eigenfunctions of (14. 2 p have the form 

M 

Hz) = lliz-a,), (5.4) 

1=1 

where {ctj \i = 1, 2, ■ ■ ■ , A/"} are roots of the polynomial which will be specified by the 
associated Bethe ansatz equations (15. 7p below. We can rewrite the Hamiltonian (14. 2 p as 

Poiz), (5.5) 



1=1 ^ ^ 



where Po{z) and Pi{z) are poljTioniials in z determined from the expansion of the products 
in fO]) . 

Dividing the Hamiltonian equation Hi/j = Ei/j by "0 gives us 

The Lh.s. of (15. 6p is a constant, while the r.h.s is a meromorphic function in z with at 
most simple poles. For them to be equal, we need to eliminate all singularities on the 
r.h.s of f l5.6p . We may achieve this by demanding that the residues of the simple poles, 
z = ai,i = 1,2, should all vanish. This leads to the Bethe ansatz equations for the 
roots : 

2^ 7 Z ^ 1 Z r + Pi(«^) = 0, 

^ = 1,2, ■■■,AA. (5.7) 

The wavefunction "0(2;) (15. 4p becomes the eigenfunction of H (14. 2p in the space V provided 
that the roots {ai} of the polynomial ^^{z) (15. 4p are the solutions of (15. 7p . 

Let us remark that the Bethe ansatz equation (15. 7p is the necessary and sufficient 
condition for the r.h.s. of (15.60 to be independent of z. This is because when (I5.7p is 
satisfied the r.h.s. of (15. 6 p is analytic everywhere in the complex plane (including points 
at infinity) and thus must be a constant by Liouville's theorem. 

To obtain the corresponding eigenvalue E, we consider the leading order expansion of 
^(^), 

^[z) = -z^-^^ai + ---. 

i=l 

It can be directly shown that the V±flil){z) have the expansions 

V+ip = -z^g (f[{2j -p-i + l-r{Af -1)) 
\i=i 



VoiJ = z"^ (^Af+^-K^+---. 

Substituting these expressions into the Hamiltonian equation (15. ip and equating the z-^ 
terms, we arrive at 

1=1 \ \ ^=1 M=* / 
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+5-' {rAf - j +pY - g 

M ki 



Y[{2j -p-t + l-r{Ar -1)) 



HHk, (A + q,-Af+l 

i=l u=l ^ 



(5. 



i=l 



where {ctj} satisfy the Bethe ansatz equations fl5.7p . This gives the eigenvalue of the 
Hamiltonian (11. ip with the corresponding eigenfunction ip{z) (15.40 . 



6 Explicit examples 

In this section we give explicit results on the Bethe ansatz equations and energy eigen- 
values of the Hamiltonian (II. ip for special cases which correspond to some established 
models frequently studied in the field of atomic and molecular physics, condensed matter, 
nuclear physics and quantum optics. 

A. Two-site Bose-Hubbard model 

This model corresponds to the special case with M = 0,r = l,s = 2 and its Hamilto- 
nian takes the simple form 

H = g'Jl + g{J+ + J_). (6.1) 

This model has been widely employed in the context of Josephson-coupled Bose-Einstein 
condensates via the realization of J±fl in terms of two bosons, 7+ = h\h2, J_ = hih\, Jq = 
|(6|6i — &2^2) (see e.g. [5] and references therein). Exact solutions of the model in terms 
of algebraic Bethe ansatz methods were first studied in [T7]. From the general results in 
the preceding section, in this case we have k = 0,p = 0, and Af = 2j. Thus (16. ip takes 
the form 

where 

P2{z) = g'z\ 

P,{z) = g'z{l-2j) + g{l + z'), 

Po{z) = g'f-2jzg. 

The Bethe ansatz equations are given by 

^ 2 a,g'il-2j)+gil + al) 

and the energy eigenvalues are 



Oii- 
i=l 
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This exact solution is equivalent to described in [TBI. 



B. Lipkin-Meshkov-Glick model 

This model is the special case corresponding toM = 0,r = 2,s = l. The Hamiltonian 
is given by [6] 

H = g'Jo + g{Jl + Jl) (6.2) 

and continues to be studied extensively (see e.g. [H] and references therein). Exact 
solution via the algebraic Bethe ansatz method is discussed in [20l[2T]. Specializing the 
general results of the preceding section to this case, we have p = 0,1, k = and JV = 
2j-p-A ^^^Yi X = 0, 1 so that is a non-negative integer. The differential operator 
representation of the Hamiltonian (16.21) is thus 

where 

P2{z) = Agz^ + Agz, 

Pi{z) = g{Q + Ap-%3)z^ + 2g'z + g{2 + Ap), 

Po{z) = gz{2j -p){2j -p-1) + g'{p- j). 

The Bethe ansatz equations are given by 

1^ 2 g{3 + 2p- A)al + g'a^ + g{\ + 2p) 



Oii - (a^ + 



^ = 1'2, ■ ■ ■ , 



and the energy eigenvalues are 

E = g'{j-X)-g{X + l){X + 2) ^ a,. 



2j-p-A 
2 



i=l 



C. Molecular asymmetric rigid rotor 

Up to an additive constant, this model corresponds to the special case with M = 0, r = 
s = 2. This is shown as follows. The Hamiltonian of the rigid rotor has the following 
form in terms of the su{2) generators J^, Jy and 

H = aJl + hJl + cJl 
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where a, b, c are constants. The model has previously been discussed in P2] as a Hamil- 
tonian which is solvable by algebraic Bethe ansatz methods. The Hamiltonian can be 
rewritten as 

H= Jq + ^^{J+ + J-) + (6.3) 

where C is the Casimir element of su(2). This shows that the molecular asymmetric rigid 
rotor is indeed a special case of (11. ip . Note that the Hamiltonian (16. 3 p of the rigid rotor 
almost has the same form as that of the Lipkin-Meshkov-Glick model. To our knowledge, 
this connection has not been noted previously. 

Specializing the general results in the preceding section to this case, we have p = 0,1, 
K = and = ^''~2~'^ ) where A = 0, 1 as in the case of the Lipkin-Meshkov-Glick model. 
The differential operator representation of the Hamiltonian (16. 3p is thus 

H = P,{z)^^ + P,{z)^ + Po{z) 
az"^ az 

where 

P^{z) = {a-b)z^ + 2{2c-a-b)z'^ + {a-b)z, 

Pi{z) = ^(3 + 2p-4j>2 + 2(2c-a-6)(l+p-j> + ^(l + 2p), 

a — b, , 2c — a — b, a + b , 

Poiz) = -^{2j-p){2j~p-l)z+ ^ {p-jf + ^jh + l)- 

The Bethe ansatz equations are given by 

2 (a-&)(3 + 2p-4)a2 +4(2c-a-6)(l+p-i)a^ + (a-6)(l + 2p) 

2(a - b){al + a^) + 4(2c - a - b)al ' 



, . 2J-P-X 
= 1'2, ■ ■ ■ , 



and the energy eigenvalues are 



2 



2c — a — 6, ,,2 a + b , , a — b,^ , 
E = (j-Af + -— j(j + l)-^(A + l)(A + 2) a,. 



4 

i=l 
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D. Tavis-Cummings model 

This model corresponds to the special case when M = r = s = ki = l. The Hamilto- 
nian is given by 



H = WiNi + g'Jo + g J+ai + J_a\^ . 



This is one of the widely studied models in quantum optics and had been exactly solved 
via the algebraic Bethe ansatz approach [23l[2l]. Applying the results in the preceding 
section gives qi = 1, p = and Af = min{2K + j — l,2j}. The differential operator 
representation of the Hamiltonian is 

where 

P2{z) = 9z\ 

P,{z) = -g(3j + 2K-2)z'' + ig'-Wi)z + g, 

Poiz) = w,{2K + j-l)-g'j + 2gjzi2K + j-l). 

The Bethe ansatz equations read 

^ 2 g{3j + 2K-2)al-{g'-w)a,-g 

and the energy eigenvalues are 

E = wi{2K + j-^-l)+g' {M-j)-g{2j-M+l) (2/s: + j - AT) J] a,, 

i=l 

where M = min{2/« + j — 1, 2j}. 

E. Two-mode generalized Tavis-Cummings model 

Finally, we consider the case when M = 2, r = s = ki = k2 = I- This gives the 
Hamiltonian 



H = wiNi + W2N2 + g'Jo + g (^J^alal + J+aia2 



which belongs to the class of su{l, 1) generalized Tavis-Cummings model discussed in 
[To] . Applying the results in the preceding section gives qi = q2 = I, p = and Af = 



min{^^^2^ — 1 + j, 2j}. The differential operator representation of the Hamiltonian thus 
reads 

12 



where 



Ps{z) = -gz\ 



with 



Pi{z) = Az^ + {g' -wi-W2)z + g, 
Po{z) = zB + F 



( f 9 

A = g {-9jK + 10j + 6K + j-5f -A- -K^ 

F = + ^2) ~ 1 + + ^(^1 ~ "^2) - g'j- 

The Bethe ansatz equations assume the form 



/3 = 1,2,---,AA 



and the energy eigenvalues are 



-^(2i-AA+l) 
where = min{^^ - 1 + j, 2j}- 



3«; . /2 



i=l 



7 Discussions 

We have derived the exact solutions of a family of Hamiltonians with the following general 
form, 

H = F{Qo) + g{Q+ + Q-) (7.1) 

whereby Q±fi are particular polynomial deformations of the sl{2) Lie algebra and F{Qo) 
is a polynomial function of Qq with real coefficients. We have seen that via the differential 
operator reahzation of these algebras, the block diagonal sectors of the Hamiltonians can 
be reahzed as higher order quasi-exactly solvable differential operators. The eigenvalues 
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of the Hamiltonians in these sectors have been obtained via the functional Bethe ansatz 
approach. 

Specific cases of the general Hamiltonian have previously been solved via the algebraic 
Bethe ansatz approach [T7 1I20I - I21] as mentioned earlier. Comparing both methods, it 
appears that the functional Bethe ansatz approach has some advantage over the algebraic 
Bethe ansatz by requiring less algebraic machinery. This advantage manifests itself in the 
fact that we have been able to give a unified solution for (11. II) through (15. 7115. 8p . Such a 
unified solution presently appears beyond the limits of algebraic Bethe ansatz approaches 
which treat the models on a case-by-case basis. It would therefore be interesting to see 
whether other classes of exactly solvable models can be easily handled by the functional 
Bethe ansatz approach. 

One avenue for further work would be to generalize the functional Bethe ansatz ap- 
proach to solve for other classes of Hamiltonians, such as g-deformed versions of the 
models discussed above. It would also be worthwhile to explore the role of polynomial 
algebra structures in connections between higher order ODEs and integrable models i.e. 
the ODE/IM correspondence [25]. 

Acknowledgments: This work was supported by the Australian Research Council. 
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